Energy Relaxation of Hot Dirac Fermions in Graphene 
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We develop a theory for the energy relaxation of hot Dirac fermions in graphene. We obtain a 
generic expression for the energy relaxation rate due to electron-phonon interaction and calculate the 
power loss due to both optical and acoustic phonon emission as a function of electron temperature 
Tc and density n. We find an intrinsic power loss weakly dependent on carrier density and non- 
vanishing at the Dirac point n — 0, originating from interband electron-optical phonon scattering 
by the intrinsic electrons in the graphene valence band. We obtain the total power loss per carrier 
- 10"^^ - 10"^ W within the range of electron temperatures ~ 20 — 1000 K. We find optical 
(acoustic) phonon emission to dominate the energy loss for To > (<)200 — 300 K in the density 
range n = 10^^ - 10^^ cm~^. 

PACS numbers: 72.20.Ht, 73.63.Bd, 63.20.kd 
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Graphene is a two-dimensional (2D) plane of carbon 
atoms in a hexagonal lattice, with the motion of the elec- 
trons near the Brillouin zone K point (the 'Dirac point') 
obeying the massless Dirac equation. This unusual quasi- 
relativistic property implies a linear energy dispersion 
^kx — hXvk of the conduction band (A — 1) and the 
valence band (A = —1) intersecting at the Dirac point, 
with a constant Fermi velocity v = lO^ms"^. The re- 
cent experimental discovery of this new 2D material has 
sparked enormous interest in understanding its funda- 
mental transport properties; for low electric fields, carrier 
transport is determined by the behavior of the carriers 
close to the Fermi level, and has been quite extensively 
studied both experimentally and theoretically [l,]. 

What still remains as an interesting open arena for ex- 
ploration is the carrier transport in graphene under high 
electric fields and in non-equilibrium (i.e. "hot-electron" ) 
conditions. With the application of a high electric field, 
carriers gain energy at a rate much faster than that for 
the carriers to lose energy to the lattice, creating a non- 
equilbrium carrier population which subsequently comes 
to an internal thermal equilibrium (among the carriers 
themselves but not with the lattice) through carrier- 
carrier scattering. These carriers (called hot carriers), 
having a higher temperature Te than the lattice temper- 
ature Tl, will then relax towards a thermal equilibrium 
with the lattice by losing energy to the lattice. As trans- 
port properties at high fields are determined by these hot 
carriers, a quantitative understanding of this hot carrier 
energy relaxation process is therefore of central impor- 
tance in high-field carrier transport; furthermore, it is a 
crucial issue affecting the performance characteristics of 
ultrafast, high-field devices. 

Experimentally, this energy relaxation process can be 
studied, following excitation with an ultrashort laser 
pulse, by measuring the transmission spectrum 0, |3, Q 
(pump-probe spectroscopy) or the luminescence spec- 



trum [H] (photoluminescence spectroscopy) from the sam- 
ple. In particular, the luminescence spectrum yields di- 
rectly the temperature of the hot carriers, and when com- 
bined with electrical (instead of optical) excitation of the 
carriers with a steady-state electric current, can provide 
a direct means to measure the carrier temperature as a 
function of the energy- loss rate @ . 

In this Letter, we develop a theory for the energy re- 
laxation of hot electrons in graphene, incorporating the 
energy loss channels due to optical phonon emission and 
acoustic phonon emission. Note that electron-electron 
interaction can only lead to energy equilibration within 
the carrier system (i.e. establishing an electron temper- 
ature), but cannot contribute to the overall energy loss 
from the carriers to the lattice. We derive an expres- 
sion for the energy relaxation rate (i.e., power loss) of 
hot electrons in graphene due to electron-phonon (e-ph) 
interaction, and obtain the power loss as a function of 
electron temperature and density. We find an intrinsic 
power loss that is present at all doping levels and does 
not vanish at zero doping arising entirely due to inter- 
band electron-optical phonon scattering by the intrinsic 
electrons in the valence band; and an extrinsic power 
loss which scales with doping density. We also find that 
the power loss is predominantly due to acoustic phonon 
emission below ~ 200 — 300 K (depending on the doping 
density) with optical phonon emission taking over as the 
dominant energy loss mechanism above this temperature 
range. The temperature throughout this paper refers to 
the electron temperature with the lattice assumed to be 
held at a low temperature. 

The rates of change of the electron and hole distri- 
bution functions describing the electron-phonon scatter- 
ing are given by (r and t dependence of the distribution 
functions are suppressed for clarity, and we set h = 1 
throughout unless specified): 
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where /i, A = ±1 is the chhahty index for electron 
(/i, A = 1) in the conduction band or hole (/i, A = —1) 
in the valence band, fk^ is the electron — 1) or hole 
(/i = —1) distribution function 7], S,k\ — ^vk — fi^ is the 
quasiparticle energy rendered from the chemical potential 
He Nq is the phonon distribution function at the phonon 
energy usq, a'^\_^ is the e-ph coupling strength which, for 
graphene, has a non-trivial chiral and momentum depen- 
dence arising from the graphene band structure and the 
e-ph interaction vertex In Eq. the first term 

within the brackets on the right hand side describes the 
scattering of an electron in the state {k,fi) into another 
state (fc — q, A) via the emission (i.e., — ^k-qX = ^^q) 
or absorption {S,k-q\ — CfeM = ^q) of ^ phonon, and the 
second term follows with an analogous physical meaning. 
The sum over A = ±1 takes into account both intra- 
band (A = /i) and interband (A ^ fi) scattering processes 
through phonon emission or absorption. 

We now state the assumptions for our model: (1) Af- 
ter the initial rapid carrier-carrier scattering, the elec- 
tron gas has established an internal thermal equilibrium 
at an electron temperature Tg described by the Fermi 
distribution function fk^ = «F(^fe/^), where npie) — 
l/[exp(e/fcBTe) -I- 1] is the Fermi function. (2) We take 
into account the fact that the emitted optical phonons 
can decay into low-energy acoustic phonons due to an- 
harmonic phonon-phonon scattering, which is character- 
ized phenomenologically by a finite optical phonon life- 
time Tah. (3) Acoustic phonons emitted by the elec- 
trons or produced through the decay of optical phonons 
thermalize immediately with the lattice (acting as a 
heat bath) which is maintained at a lattice temperature 
Tl < Te- Recent ultrafast optical spectroscopy experi- 
ments on graphene [1, |3| find that the time for the elec- 
trons to equilibrate among themselves spans ~ 100 fs, 
and the subsequent thermalization of the electron gas 
with the lattice lasts for ~ 1 ps. It is within this pi- 
cosecond time scale that the electrons lose most of their 
energy through e-ph scattering as electron-electron scat- 
tering does not dissipate energy from the electron gas as 
a whole. 

We first ignore the effect of a finite phonon lifetime by 
taking Tah — 0, assuming that the emitted phonons im- 
mediately thermalize with the lattice and all the energy 
lost from the electrons to the phonons is also immedi- 
ately lost to the lattice. The energy loss rate dEkf^/dt 
of a single electron with momentum k can be obtained 
from Eq. ([T]) by inserting the energy change S,k-q\ — ^kfi 



of the electron due to scattering under the q integral. 
The total energy loss rate P of the entire system of 
electrons then follows by summing the resulting single- 
electron energy loss rate over all states of momentum k 
and chirality /i, taking into account the degeneracy fac- 
tors due to spins gg = 2 and valleys = 2 to give P = 
~5vffs J2fj. 12k d-Bfc^/di. Awaking use of the integral iden- 
tities d{Ck,i-^k-q\T^q) = J dujS{^k^-S,k-q\+(^)S{ujq±Uj) 
in the expression of dEkp,/dt obtained above, we arrive, 
after some algebra, at the following expression for the 
total power loss from the electrons: 
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where (lu) — l/[exp(cL;/fcBrc^L) — 1] stand for the 
Bose distribution functions evaluated at the electron 
Tc and lattice Tl temperatures respectively, T>{q, uj) = 
2ujq/ {uj'^ — ujq + iO'^) is the phonon Green function. 
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is the phonon self-energy [lOI| at the electron temperature 
Tc, and 'Im' in Eq. ([2]) stands for the imaginary part. 
Eq. © generalizes the Kogan formula [Il[ for the power 
loss in an e-ph coupled system widely used in regular 
metals and semiconductors [12| with a parabolic energy 
band to a chiral two-band system (to which graphene 
belongs as a special case), embodying both intraband and 
interband electronic transitions as well as the non-trivial 
chiral and momentum dependence of the e-ph coupling. 

We first consider the energy relaxation due to opti- 
cal phonon emission. The LO phonon mode at the Bril- 
louin zone center F in graphene is characterized by the 
phonon energy ujq = 196 meV and e-ph coupling ct'^\._q = 
glp[l - nXcos{(l3k + (kk-q - 2^]/2 [i,[i,[l3|, where g^^ is 
the e-ph couphng constant [l3( and (jjk = ta,ii^^ {ky / kx) ■ 
The imaginary part of the phonon self-energy Eq. ([3]) de- 
scribes the damping of the phonon mode due to electron- 
hole pair excitations. For the F point optical phonons 
in graphene, nP*^ is different from the graphene polariz- 
ability ^14] 11 due to the different chiral and momentum 
dependence in the e-ph coupling. We have obtained an 
exact ana,lytical expression for the ImnP*^ at zero tem- 
perature [l5|; the finite-temperature phonon self-energy 



is then obtained from the zero-temperature expression as 
0: 
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here we have written out the dependence of ImllP'^ (i.e. 
the imaginary part of Eq. ([3])) on the temperature T and 
the chemical potential Pc{T) explicitly for clarity. Pc{T) 
is determined by requiring that the integral over all the 
electronic states of the Fermi function gives the electron 
density n = /dez^(£)/{exp[(e — /Xc(re))/fcBT'e]-|-l}, where 
v{e) = g^ggS /2t:v'^ is the graphene electronic density of 
states. Unlike in regular 2D parabolic system, a closed 
form for p, cannot be obtained in the case for graphene, 
and one must numerically solve for the root of the equa- 
tion which results from the integration of the Fermi func- 
tion: n = -[2(fcBTc)^/7rw^]Li2[— exp(/ic/A;Bre)], where 
Li2{z) — X]fcLi(^'^/^^) is the dilogarithm function. 

In Eq. ([3]), we can identify two contributions to the 
phonon self-energy at zero temperature IIPp = IlPP-t-IlPP; 
with nPP originating from the extrinsic carriers and 
therefore dependent on the Fermi level; IIPp from the 
intrinsic electrons in the valence band (i.e., the 'Dirac 
sea') and independent of the Fermi level. For optical 
phonons, both the extrinsic and the intrinsic parts con- 
tribute to the total power loss P = Pcxt + -Pint, and as 
the extrinsic carrier density n is tuned to zero, Pext goes 
to zero, but Pint remains finite even at zero extrinsic car- 
rier density n = 0. The intrinsic part of the phonon 
self-energy ImllPP has a simple form, which allows for an 
exact analytic derviation of the intrinsic power loss Pint 
from Eqs. ©-(gl): 
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here N = nA {A is the sample area) is the number of 
extrinsic carriers and (ffop)^ = ffop^/^^'*'^ the dimension- 
less e-ph coupling P. Physically, Pint correponds to the 
power loss due to optical phonon emission through in- 
terband transitions of the valence band electrons, and is 
only weakly dependent on the extrinsic carrier density 
through pc in the Fermi function. We find that this in- 
trinsic power loss is not small, and for To = 300 — 700 K 
at rt = lO^'^cm"^, the intrinsic power loss per unit area 
PintM ~ 10^ - 10^ Wm-2. _ 

Recent experiments ^17] and theory [iS] show that 
Tah for graphene is of the order of picoseconds, and 
therefore anharmonic phonon-phonon scattering occurs 
at a comparable (and slower) rate than e-ph scatter- 
ing, causing an accumulation of non-equilibrium opti- 
cal phonons (known as 'hot phonons'). Some of these 
hot phonons are then reabsorbed back by the elec- 
tron gas, thus reducing the overall power loss. We 
now take this effect into account by incorporating [l^ 




0.001 0.003 0.005 0.007 0.009 0.011 

1/T (k') 

FIG. 1: (Color online) Power loss per carrier versus inverse 
electron temperature l/Te for optical phonons (black) and 
acoustic phonons (red) at different electron densities n = 10^^ , 
10^^, 10^^ cm~^. The slopes of the curves for optical phonons 
are approximately given by luq at low To. 



a finite phenomenological Tah in our theory, and ob- 
tain the power loss as P = ^TlqS dw(u;/7r)[n^(w) — 
n^(w)]lml?((7, w)/[r(q, oj) + Tah], where we have written 
= — 2ImnP^((7, Ld) as the phonon damping rate due to 
electron-hole pair excitation. At lower electron tempera- 
tures, electrons do not have enough energy to emit high 
energy optical phonons, and acoustic phonon emission 
becomes the dominant energy loss mechanism. The LA 
mode acoustic phonon in graphene at the F point has an 
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energy dispersion ujq — CgQ and e-ph coupling [8| a 
gl,{q)[l + pXcos{(j)k^q - (/)fc)]/2, wher e Cs = 2 x 10'' ms"^ 
is the phonon velocity, (?ac(9) = Dq^h/ pAujq is the cou- 
pling constant with p being the graphene mass density, 
and D — 16 eV the deformation potential. In contrast 
to the case of optical phonons, we note that a'l\_q for 
acoustic phonons has the same chirality and momentum 
dependence as in the graphene polarizability, and there- 
fore the phonon self-energy for LA phonons is given from 
Eq. © simply by nPP((7,u;) = 5^(0)0(9, w), withn(q,tj) 



being the polarizability in Ref. [14 1 



As in the case for optical phonons, the acoustic phonon 
self-energy contains an extrinsic contribution IIPp and an 
intrinsic contribution IIPP; however, we find that the in- 
trinisc part does not contribute to the power loss and 
the power loss due to acoustic phonon emission origi- 
nates entirely from the extrinsic contribution P = Pcxt- 
Physically, the intrinsic contribution corresponds to in- 
terband electron transitions across the conduction and 
the valence bands. The acoustic phonon mode, having 
an energy ojq = Cgq with Cg smaller than the graphene 
band velocity v, does not provide a possible channel for 
interband transition, which requires an energy greater 
than vq. Emission of acoustic phonons is therefore only 
possible through intraband transitions. 

We now calculate the power loss Popticai due to optical 



phonons (we use r^h = 3.5 ps from Ref. [l^) and Pacoustic 
due to acoustic phonons as a function of Te and n, with 
IJLc{Te) determined at each Te and the finite-temperature 
phonon self-energy obtained by evaluating Eq. The 
lattice temperature is taken as zero Tl — 0. Fig.[l]shows 
the power loss per carrier P/N from the two contribu- 
tions versus inverse electron temperature I/Tq. The ap- 
proximate exponential behavior of Popticai versus l/Te is 
due to the fact that the electrons capable of emitting 
optical phonons have an amount of energy higher than 
Wo, which lie in the high-energy tail of the Fermi dis- 
tribution with a population ~ exp(— wo/fcBTc)- As the 
phonon energy is quite high ojq = 196 meV in graphene, 
the power loss through optical phonons decrease with Tq 
about an order-of-magnitude faster than that in GaAs 
(where ujq — 36meV). On the other hand, Pacoustic de- 
creases with Tc much more slowly, lying within the range 
of 10-12 - 10-9 W for To = 100 - 1000 K. 

The total power loss is given by the sum of the contri- 
butions from the optical phonons and acoustic phonons 
(Fig. [2|). At small values of 1/To (high temperatures), 
the energy loss is predominantly through optical phonon 
emission, with the power loss behaving approximately 
exponentially. Topticai decreases as temperature is de- 
creased, and the power loss through acoustic phonons 
-Pacoustic becomes increasingly important. The crossover 
of the energy loss from predominantly optical phonon 
emission to predominantly acoustic phonon emission de- 
pends on the electron density, occurring at an increasing 
temperature with density To ~ 200 K for n — 10^^ cm-^, 
- 250 K for 10^2 cm-2 and - BOOK for IQi^cm-^. 
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FIG. 2: (Color online) Inverse electron temperature 1/Tc ver- 
sus the total power loss per carrier P/N. 1/Tc is plotted 
against P/N since P is the input variable in experiment (the 
power loss P to the lattice must be equal to the experimental 
power input to the electrons) while Te is the observed out- 
put. Inset: Total power loss per carrier shown within a wider 
temperature range down to Te = 20 K. 

We show in Fig. [3] the total power loss as a func- 
tion of electron density n for different values of To. At 
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FIG. 3: (Color online) Total power loss per carrier P/N versus 
electron density n at different values of Tc. 



To ~ 100 K, the energy loss is mainly due to acoustic 
phonon emission, and P/N increases with density. For 
Tc ^ BOOK, P/N shows an upturn as n is reduced to- 
wards zero, reflecting the portion of intrinsic power loss 
coming from optical phonon emission which must give 
P/N^(xa.sn^O since Pi„t is finite at ri = 0. 

In conclusion, we emphasize that the intrinsic power 
loss has its origin from the presence of the Dirac sea 
in graphene, and therefore energy relaxation resulting 
in such intrinsic power loss occurs not just for undoped 
graphene with n — 0, but for doped graphene at all car- 
rier densities n > as well. We find that the total power 
loss per carrier taking account of both optical and acous- 
tic phonon emission lO^^^ — 10"^ W for electron tem- 
peratures To 20 — 1000 K. Our results obtained for 
Tl = should remain valid as long as Tl ^ To is satis- 
fied in experimental situations. 
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